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Abstract: We construct a family of pairwise commuting operators such that 
the Jack symmetric functions of infinitely many variables Xi,x%, ... are their 
eigenf unctions. These operators are defined as limits at N — > oo of renormalised 
Sekiguchi-Debiard operators acting on symmetric polynomials in the variables 
x±, . . . ,Xn- They are differential operators in terms of the power sum variables 
p n = Xi + X2 + ■ ■ ■ and we compute their symbols by using the Jack reproducing 
kernel. Our result yields a hierarchy of commuting Hamiltonians for the quantum 
Calogero-Sutherland model with infinite number of bosonic particles in terms of 
the collective variables of the model. Our result also yields explicit shift operators 
for the Jack symmetric functions. 
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1. Introduction 

1.1. Calogero- Sutherland model. This quantum model describes a system of N 
bosonic particles on a circle IR/7rZ with the Hamiltonian [5,19,20] 

d 2 , ^ W-i) 



9qf fZj sin fe-^) 

where ^ qi, • • • ,qN < ^- Being translationally invariant Hqs commutes with 
the momentum operator 



After eliminating the vacuum factor 



and then passing to the exponential variables Xj = exp(2i(/j) and the parameter 
a = more common in the mathematical literature, the Hamiltonian becomes 

-1 rr 2 tt(2 ) N 3 ~N 

u oHcqou = — H^' + 



where 



a JV 6^ 



Respectively, Pes gives rise to the operator 



The operators and commute and act on symmetric polynomials of 
the variables x\, . . . , xn ■ It is known that both operators can be included into 
a quantum integrable hierarchy, that is into a polynomial ring of commuting 
differential operators with N generators of orders 1, . . . , N respectively, called the 
Sekiguchi-Debiard operators [6,8,16]. The Jack symmetric polynomials [8] are 
joint eigenfunctions of the hierarchy. They are labelled by partitions of 0, 1, 2, . . . 
with no more than N non-zero parts. 

In combinatorics it is quite common to extend various symmetric polynomials 
to an infinite countable set of variables. These extensions are called symmetric 
functions. In particular, the extensions of the Jack symmetric polynomials are 
well studied [8]. They are labelled by partitions of 0, 1, 2, ... . However, no explicit 
expressions for higher commuting Hamiltonians corresponding to the infinite set 
of variables have been yet available in the Jack case, with an exception of a few 
lower order operators. The main purpose of the present article is to fill up the gap, 
by studying the limits of the Sekiguchi-Debiard operators as iV -> oo and giving 
explicit expressions for the resulting commuting Hamiltonians. 
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As another application of our result, we construct elementary shift operators 
for the Jack symmetric functions. In terms of the labels, our operators correspond 
to decreasing any given non-zero part of a partion by 1 and to the operation on 
partitions inverse to that, see our formulas (2.27) and (2.25) respectively. For the 
origins of this construction see the work [18] and references therein. 



1.2. Collective variables. The standard way to treat the N = oo case is to rewrite 
the Hamiltonian (1.1) in terms of the power sums (the "collective variables" in 
the condensed matter physics terminology) 

Xi+.-. + Xjy where n = l,...,N (1.2) 

and to take the limit at N — > oo afterwards. Denoting the limit of the power 
sum (1.2) by p n where n = 1,2,... the resulting Hamiltonian reads [9] 



H^= lim a(H^-NH^) 



(2) _ Y// !ii , 

c/ 2 



S~] (a(m + n)p m p n —^- \-a 2 mnp m+n 



n=l 



dpmdpr, 

QPr, 



+ 



°° d 
(a-1) ^2 an 2 p n - — . (1.3) 



The first and the second summands in the middle line of the above display are 
known as splitting terms and joining terms respectively. 

Consider the vector space A = C[pi,p2> • • • ] and equip it with the operators 
a n where n G Z \ {0} , defined on the polynomials / G A by 



«n(/) = 



P-nf if n<0, 
andf/dp n if n>0; 



see also the equality (2.14) below. The operators a n satisfy the relations 

[a m ,a n ] = ma5 m+n fi ■ (1.4) 

Thus A becomes a highest weight module for an infinite-dimensional Heisenberg 
Lie algebra. In terms of the operators a n the Hamiltonian (1.3) takes the form 

oo oo 

H {2) _ ^ (a- m a- n a m+n + a- m - n a m a n ) + (a-l)^2na- n a n . 

m,n=l n=l 

Similarly, the yields a first order differential operator commuting with H ^ 

OO q oo 

= lim aH$ = anp n - — = a_ n a n . (1.5) 

n=l Ln n=l 
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1.3. Quantum field theory. In terms of the field cp(s) on the circle S 1 = M/2nZ 
(p(s) = a n exp( — ins) for s <E S 1 

the commutation relations (1.4) read 

[(p(s),ip(t)} = 2nia5'(s -t) . 

Setting an = and using the Wick normal ordering : : with the operators a n 
for n < to the left and for n > to the right, the Hamiltonian becomes 



1 1 f ds 2 , 

2~a ^; W= 2^ Js^ :iP{S): 



while becomes the Hamiltonian of quantized Benjamin-Ono equation [1, 13] 

1 v a — 1 x -v 

- : aia m a n : + — — ^ 



l-\-m-\-n—0 ~ m+n=0 



where "H stands for the Hilbert transform 



(U(f)(s) =p.v. / ^ cot^y^ y?(t) 



In the particular case a = 1 the Jack symmetric polynomials degenerate into 
Schur polynomials. The Benjamin-Ono equation respectively degenerates into 
the dispersionless KdV (also called Burgers) equation. An explicit construction 
of a countable set of commuting Hamiltonians for the quantum dispersionless 
KdV can be obtained via boson-fermion correspondence and is available in terms 
of recurrence relations [12] or a generating function [14]. 

The higher quantum Hamiltonians for any parameter a are constructed in the 
present article, see the theorem in Subsection 2.5. Note that in the case a = 1 
our Hamiltonians are different from those obtained in [12,14] being rather 
their polynomial combinations, see for example the relations (2.16) and (2.17). 

In their turn, Jack symmetric polynomials can be regarded as degenerations 
of the Macdonald polynomials in the variables x\, . . . , xn also depending on two 
formal parameters q and t . The Jack case corresponds to q = t a where t — > 1 . 
The Sekiguchi-Debiard operators can be then regarded as degenerations of the 
Macdonald operators [8] acting on the symmetric polynomials in x\, . . . , xn ■ Our 
theorem generalizes to the Macdonald case [10], see also the earlier works [2, 17]. 
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1.4- Plan of the article. In the next section we recall some basic facts from the 
theory of symmetric functions and set up the notation. We then introduce the 
Jack polynomials and Sekiguchi-Debiard differential operators. Our main tool 
is the notion of the symbol of an operator relative to the reproducing kernel for 
Jack polynomials. After establishing the basics, we state our main result which 
is an explicit formula for the symbol of the generating function of commuting 
Hamiltonians. Then we explicitly construct our elementary shift operators for 
the Jack symmetric functions. We finish Section 2 with reducing the proof of our 
theorem to certain determinantal identities which are then proved in Section 3. 

In this article we generally keep to the notation of the book [8] for symmetric 
functions. When using the results from [8] we will simply indicate their numbers 
within the book. For example, the statement (1.11) from Chapter I of the book 
will be referred to as [1. 1.11] assuming it is from [8]. We do not number our own 
lemmas, propositions, theorems or corollaries because we have only one of each. 

2. Symmetric functions 

2.1. Monomial functions. Fix any field F. For any positive integer N denote 
by An the F-algebra of symmetric polynomials in N variables xi, . . . , xn . The 
algebra A N is graded by the polynomial degree. The substitution xn = defines 
a homomorphism An — > A/v-i preserving the degree. Here Aq = F. The inverse 
limit of the sequence 

A!<-A 2 <-... 

in the category of graded algebras is denoted by A . The elements of A are called 
symmetric functions. Following [8] we will introduce some standard bases of A . 

Let A = ( Ai, A2, • • • ) be any partition of 0, 1, 2, ... . The number of non-zero 
parts is called the length of A and is denoted by £{X) . If £(X) ^ N then the sum 
of all distinct monomials obtained by permuting the N variables 
is denoted by m\(xi, . . . , xn) • The symmetric polynomials m\(xi, . . . ,xn) with 
£(X) ^ N form a basis of the vector space An ■ By definition, for £(X) ^ N 

m x (x u ...,x N )= £ E ^a" 1 ^---^ (2-1) 

where we write k instead of £(X) . Here &k is the symmetric group permuting 
the numbers 1, . . . , k and 

d x = k 1 \k 2 \ ... (2.2) 
if fei, &2, . . . are the respective multiplicites of the parts 1, 2, ... of A . Further, 

/ n x j m x (xi,...,x N -i) if ^(A)<iV; 

m x {x 1 ,...,x N -i,0) = I Q . f £ ^ =N (2-3) 

Hence for any fixed partition A the sequence of polynomials m\(xi, . . . , xn) with 
iV ^ £(X) has a limit in A . This limit is called the monomial symmetric function 
corresponding to A . Simply omitting the variables, we will denote the limit by 
m\ . With A ranging over all partitions of 0, 1, 2 . . . the symmetric functions m\ 
form a basis of the vector space A . Note that if £(X) = then we set ni\ = 1 . 
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2.2. Power sums. For each n = 1, 2, . . . denote p n (xi, . . . ,xn) = + • • • + ^at • 
When the index n is fixed the sequence of symmetric polynomials p n (xi, . . . ,xn) 
with N = 1, 2, . . . has a limit in yl , called the power sum symmetric function of 
degree n . We denote the limit by p n . More generally, for any partition A put 

P\=P\ 1 P\ 2 --- (2.4) 

where we set po = 1 . The elements p> form another basis of A . In other words, 
the elements pi,j>2, • • • are free generators of the commutative algebra A over F . 

The basis of p\ can be related to the basis of monomial symmetric functions as 
follows. For any two partitions A and fi denote by the number of mappings 
9 : {1,.. .,£(n)} ->• {1,2, ...} such that 

/Uj = Aj for each z = 1, 2, . . . . (2.5) 

For any such 9 the partition fi in (2.5) is called a refinement of A . Note that if 
-Raa* 7^ then A and \i are partitions of the same number. Moreover, then by 
[1.6.10] we have fi ^ A in the natural partial ordering of partitions: 

Ati^Ai, hi + /j 2 <: Ai + A 2 , ... • 

By [1.6.9] we have 

Pv = ^R\^m\. (2.6) 



2.3. Jack functions. Now let F be the field Q(a) where a is another variable. 
Define a bilinear form ( , ) on the vector space A by setting for any A and fx 



where 



z x = l fel /ci!2 fe2 /c 2 ! 



(2.7) 



in the notation (2.2). This form is obviously symmetric and non-degenerate. By 
[Ex. VI. 4. 2] there exists a unique family of elements P\ e A such that 

(P A ,P M ) = for \^fi 

and such that any P\ equals m\ plus a linear combination of the elements m M 
with < A in the natural partial ordering. The elements P\ G A are called the 
Jack symmetric functions. Alternatively, they can be defined as follows. 
Denote by A(x±, . . . , xn) the Vandermonde polynomial of iV variables 



det 



N-j 



N 
i,J-- 



i = n ( x * 



Put 



Sn{u) = A(xi, . . . , xn) 1 • det x™ J : (u + j — 1 — axidi) 



N 



(2.8) 
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where u is a variable and di is the operator of partial derivation relative to Xi . 
Here the determinant is defined as the alternated sum 

N 

E II {xf-^iu + aW-l-axidi)) (2.9) 

where as usual (— 1) CT denotes the sign of permutation a . In every product over 
% = 1,...,N appearing in (2.9) the operator factors pairwise commute, hence 
their ordering does not matter. Further, Sn(u) is a polynomial in the variable u 
with pairwise commuting operator coefficients preserving the space An , see for 
instance [Ex. VI. 3.1]. We will call the restrictions of these coefficients to the space 
An the Sekiguchi-Debiard operators. By [Ex. VI. 4. 2] the latter operators have a 
common eigenbasis in An parametrized by partitions A of length £(X) ^ N. The 
eigenvectors are called the Jack symmetric polynomials. 

For each A with £(X) ^ N there is an eigenvector denoted by P\(xi, . . . , xn) 
which is equal to m\(xi, . . . ,xn) plus a linear combination of the polynomials 
m^xi, . . . ,xn) with \x < A and £(fi) ^ N . It turns out that each coefficient in 
this linear combination does not depend on N . Note that if A and \i are any two 
partitions of the same number such that A ^ \x , then £(X) ^ £([i) due to [1. 1.11]. 
It follows that the polynomials P\(xi, . . . ,xn) enjoy the same stability property 
as the polynomials m\(xi, . . . ,xn) in (2.3): 

p , . (P x (x 1 ,...,x N _ 1 ) if £(\)<N; 

In particular, the sequence of polynomials P\(xi, . . . ,xn) with ^ £(X) has a 
limit in A . This is exactly the Jack symmetric function P\ . The eigenvalues of 
Sekiguchi-Debiard operators acting on A N are also known. By [Ex. VI. 4. 2] 

N 

S N {u) P\(x!, . . .,x N ) = Y[(u + i - 1 - a Xi) ■ P\(x!, . . .,x N ) . (2.11) 



2.4- Reproducing kernel. In this subsection we will regard the elements of A as 
infinite sums of finite products of the variables x±,X2, ■ ■ ■ ■ For instance, we have 

Pn = x l + x 2 + • • • 

for any n ^ 1 . When we need to distinguish x±, X2, ■ ■ ■ from any other variables, 
we will write f(xi,X2,...) instead of any / e A . Now let y\, 1/2, ■ ■ ■ be variables 
independent of xi,X2,... ■ According to [VI. 10.4] with the bilinear form (2.7) 
one associates the reproducing kernel 

00 

n = Y[(i-x iyj )-^ a . (2.12) 

This II should be regarded as an infinite sum of monomials in x\,X2, ■ ■ ■ and in 
yi,V2, • • • by expanding the factor corresponding to i,j as a series at Xiyj — > . 

The property of 77 most useful for us can be stated as the following lemma. 
For any / G A denote by / * the operator on A adjoint to the multiplication by 
/ relative to the bilinear form (2.7). Note that here / = f(xi, X2, ■ ■ •) • 
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Lemma. We have 

r(II)/n = f(y u y 2 ,...). (2.13) 

Proof. The commutative algebra A is generated by the elements p n with n > 1 . 
Therefore it suffices to prove (2.13) for f = p n only. Consider the operator d/dp n 
of derivation in A relative to p n = p n (xi, £2, • • •) • By the definition (2.7) we have 

p*=and/dp n (2.14) 

On the other hand, by taking the logarithm of (2.12) and then exponentiating, 

00 

11 = eX P ( Pn{xi,X 2 , ■ ■ ■) p n {yi,V2, ■ ■ -)/any 

n=l 

The relation (2.13) for / = p n follows from the last two displayed equalities. □ 



2.5. Main result. Let F = Q(a) as in the previous two subsections. For N ^ 1 
let pn he the homomorphism An — > An-i defined by setting xn = , as in the 
beginning of Subsection 2.1. Denote 

A N (u) = S n (u)/(u)n (2.15) 

where we employ the Pochhammer symbol 

(u) N = u(u + l)...(u + N -1). 

The right hand side of the equation (2.15) is regarded as a rational function of 
u with the values being operators acting on the space An ■ Due to the stability 
property (2.10) of Jack symmetric polynomials, the equation (2.11) implies that 

p N A n (u) = A N -i(u) p N 

where Aq{u) = 1 . So the sequence of An(u) with N ^ 1 has a limit at iV — > 00. 
This limit can be written as a series 

A(u) = 1 + A^/(u) 1 + A^/(u) 2 + ... 

where the coefficients A^\A^ 2 \ ... are certain linear operators acting on A . By 
definition, the Jack symmetric functions are joint eigenvectors of these operators. 
In particular, the operators A^\ A^ 2 \ ... pairwise commute, and are self-adjoint 
relative to the bilinear form (2.7). We call them the Sekiguchi-Debiard operators 
at infinity. Due to the property (2.10) their definition immediately implies that 

A^P x = if £(\)<k. 
It is also transparent from (2.11) that for any homogeneous / G A 

A^f = -adegf. 

Hence in the notation (1.5) 

-A^=H^. (2.16) 
The operator A^ is well studied [Ex. VI. 4. 3]. In particular, it is known that 

A^(A^ + 1)-2A^ =H<V (2.17) 
in the notation (1.3). The main result of our article is the more general 



Sekiguchi-Debiard operators 



9 



Theorem. In the notation (2.2) for each = 1,2,... we have 

A (k) = { _ 1) k J- d x m x m* x (2.18) 
e(\)=k 

where X ranges over all partitions of length k . 

By inverting the relation (2.6) any monomial symmetric function m x can be 
expressed as a linear combination of the functions p^ where A , \x are partitions of 
the same number and A ^ . By substituting into (2.18) and using (2. 4), (2. 14) 
one can write each operator A^ k > in terms of p n and d/dp n where n = 1,2,... . 
In particular, one recovers the above formulas for the operators A^ 1 ' and A^ 2 \ 



2.6. Shift operators. In this subsection we will get a corollary to our theorem by 
using the following particular case of the Pieri rule for Jack symmetric functions. 
By [VI. 6. 24] for any partition \x the product piP^ equals the linear combination 
of the symmetric functions P\ with the coefficients 

n~ ! ajXi-X^-i+j-l yr a(X % - Xj - 1) - i +j + 1 
; "(A, Xj 1) i + j '11 aiK-X^-i + j ■ } 

where A ranges over all partitions such that the sequence Ai, A2, ... is obtained 
from jUi, ■ ■ ■ by increasing one of its terms by 1 and i is the index of the term. 

Further, by [VI. 6. 19] the above stated equality implies that for any partition 
A the symmetric function dP\jdp\ = a~ x p\P\ equals the linear combination of 
the P M with the coefficients 

'fr 1 o(A t -j-l) + A--z + l ^ ajXi-j + V + X'j-i 

/J - j) + Xj - i 1 = 1 a (X i -j) + \'.-i + l [ - } 

where fx ranges over all partitions such that the sequence fj,\, fj,2, ... is obtained 
from Ai, A2, • • • by decreasing one of its terms by 1 and i is the index of the term. 
As usual, here A' = (A{, A£, . . .) is the partition conjugate to A. 
Now define 

B(«) = B«/(tt) 1 + B( 2 V(«) 2 + ... 
where B^ 2 \ ... are the linear operators acting on A determined by setting 

S<*+D = (_!)* d X uim Xul m* x 
e(\)=k 

where AUl denotes the partition obtained from A by appending one extra part 1. 
Also define 

c( w ) = c«/(w) 1 + c( 2 7(«) 2 + ... 

where 

C (k+i) = { _ 1)k a -i J2 d xu im x m xul . 

e(x)=k 

Note that B(u)* = aC(u). Moreover, we will have a corollary to our theorem. 
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By the definition (2.15) of the series A(u), for any partition A we have 

^)p> = n " +i : 1 7 Ai -^. < 2 - 2i > 

u + % — 1 

1=1 

see (2.11). In the infinite product displayed above the only factors different from 
1 are those corresponding to the indices i = 1, . . . , £(X) . For any such i consider 
the product 

1 TT " + J-1-"A 3 - . (2 . 22) 
it + i - 1 1 ^ u + j - 1 v ' 

Corollary, (i) For any given partition [i we have the equality 

B{u)P lx = Y J BxM p * 

A 

where B\^{u) is the product of (2.19) by (2.22) while X ranges over all partitions 
such that the sequence Ai, A2, . . . is obtained from ^2, • • • by increasing one of 
its terms by 1 and i is the index of the term. 

(ii) For any given partition A we have the equality 

where C^\(u) is the product of (2.20) by (2.22) while [i ranges over all partitions 
such that the sequence ^2, ■ ■ ■ is obtained from Ai, A2, . . . by decreasing one of 
its terms by 1 and i is the index of the term. 

Proof. By [Ex. 1.5.3] our definition of C(u) implies that the operator commutator 

[A(u),d/dp 1 ]=aC(u). (2.23) 

Part (ii) of the corollary follows from our theorem by using the relation (2.23). 
By taking the relation adjoint to (2.23) relative to the bilinear form (2.7) we get 

{ Pl ,A(u)]=aB(u). (2.24) 

Part (i) of the corollary follows from our theorem by using the latter relation. □ 

The relations (2.23) and (2.24) were used in [15, Sec.l] as definitions of the 
series C(u) and B(u) respectively. Our theorem provides explicit formulas for 
the series defined in that way. Hence we can explicitly construct the elementary 
shift operators for Jack symmetric functions, see the equalities (2.25) and (2.27). 

Let the partition A be fixed. Then for i = 1, . . . , £(X) the elements a\i — i + 1 
of the field Q(ct) are pairwise distinct. Therefore by the part (i) of the corollary 
for the partition /i corresponding to any of these indices i we have 

B{aX t - i + 1) P M = B\^(aXi -i + l)P x (2.25) 
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where the coefficient B\^(a\i — i + 1) is the product of (2.19) by 

£(\) £{X) 

U aX ._ l+1 - Uiali-aXj-i + j). (2-26) 
j=1 aA t i + j j=1 

The left hand side of the equality (2.25) should be understood as the value in A 
of the rational function B(u) at the point u = a\i — i + 1. Similarly, by (ii) 

C(aX l -i + l)P x = C^ x (a\ l - i + 1) P M (2.27) 

where the coefficient C^\(a\i — % + 1) is the product of (2.20) by (2.26). 



2. 7. Reduction of the proof. In this subsection we reduce the proof of our theorem 
to proving a certain determinantal identity for each N = 1,2,... . This identity 
will be proved in the next section by using the induction on N . 

By the lemma from Subsection 2.4 our theorem is equivalent to the equality 

°° (-l) k 

A{u){n)/n =Y, W- d x mx(x 1 ,x 2 ,...)mx(yi,y2,...) (2.28) 

k=0 ^ U * k £(\)=k 

where the coefficients of the series A(u) are regarded as operators acting on the 
symmetric functions in the variables xi, X2, ■ ■ ■ ■ Here we set (u)o = 1 . It suffices 
to prove for each N = 1, 2, . . . the restriction of the functional equality (2.28) to 

xn+i = x N+2 = . . . = 0. (2.29) 

By the very definition of A(u) the restriction of the left hand side of (2.28) 
to (2.29) as of a function in the variables xi,x 2 , ■ ■ ■ equals 

A N ( u )(n N )/n N (2.30) 

where we denote 

N oo 

n N = HI[{i-x i y s )- 1 ' a . 

i=l j=l 

Simply by the definition of the monomial symmetric function mx(xi,X2, • • •) its 
restriction to (2.29) is mx(xi, . . . , xn) if ^(A) ^ and vanishes if £(X) > N. 
Therefore the restriction of the right hand side of (2.28) to (2.29) equals 

N (-l) fc 

^2 ( \ S dxmx{x 1 ,...,x N )mx{yi,y2,---)- (2.31) 

k=o ^ u ' k e(\)=k 

Further, due to [VI. 2. 19] to prove the equality of (2.30) to (2.31) it suffices to set 

Vn+i = yN+2 = . . . = . 



However, we will keep working with the infinite collection of variables y±, y 2 , . 
This will simplify the induction argument in the next section. 
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Let us compute the function (2.30). It depends on the variable u rationally. 
It is also symmetric in either of the two collections of variables xi, . . . , xn and 
yi, y2, . . . . This function can be obtained by applying to the identity function 1 
the result of conjugating An(u) by the operator of multiplication by i7jv- 

Conjugating the operator (2.9) by IIn amounts to replacing each di in (2.9) 
with the sum 

oo 

Here we are just adding to each di the logarithmic derivative of the function IIn 
relative to Xi . Hence conjugating (2.9) by the multiplication by IIn yields 



jv oo 

e (-!)• n k- co («+'W -i -<.*.*-£ ^f-)) 

aee N i=i i=i iyi 



Here in any single summand each of the factors corresponding to i = 1 , . . . , N 
does not depend on the variables Xj with j ^ i . Therefore when applying the 
latter operator sum to the identity function 1 we can simply replace each di with 
zero. Then we get the function 



N oo 

E(-'rn(^ ) (« + ^-i-Ei^)) = 



det 



+ ( 2 .32) 



1=1 X ^ 



It follows that the function (2.30) is equal to the determinant (2.32) divided 
by the Vandermonde polynomial A(x\, . . . , Xn) and by the Pochhammer symbol 
(u)n, see (2.8) and (2.15). This ratio is equal to the right hand side of (2.31) by 
the proposition in Subsection 3.1, see the argument at the end of that subsection. 
We will prove the proposition in Subsections 3.2 to 3.4. Thus we will complete the 
proof of our theorem. Note that another proof of this theorem can be obtained 
by using the results of [2, Sec. 3] and [17, Sec. 9] on the Macdonald operators. 



3. Determinantal identities 

3.1. Getting the theorem. Let F be any field. Consider the rational function of 
two variables u, v 

V(u,v) = -^— (3.1) 
uv — 1 

with values in F . This function is a solution of the equation 

(u - v) \P(u, w) \P(u, w) =u \P(v, w) - v \P(u, w) . (3.2) 

Here w is a third variable. Any non-zero rational solution \P(u,v) of (3.2) is 

u 

u — ip(v) 
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where ifi(v) is an arbitrary rational function of a single variable. In particular, 
here by choosing ip(v) = v~ l we get the solution (3.1). 

Let xi, . . . , xn and yi, y2, . . . be independent variables. Here we assume that 
N ^ 1 . In the next three subsections we will prove the following proposition. 

Proposition. For any solution &(u, v) of the equation (3.2) we have an identity 



det 



oo 

Xl N - j (u + j-i+J2 <p(*i> yi))] l 3=1 



1=1 

N 



A(x u x N ) ^2 (u + k) . . .(u + N - 1) II *( x ^Vjr) (3-3) 

k=0 r=l 

ji, — ,3k 

where all the indices ii, . . . ,ik G {1,..., N} are distinct, the indices ji, ■ ■ • ,jk £ 
{1,2,...} are all distinct too, and the sum is taken over all collections of these 
indices such that different are all the corresponding sets of k pairs 

{(zi,ji),...,(ifc,j fc )}. (3.4) 

For any k ^ 1 take the symmetric group &k • Using the permutations a G &k 
the sum over the indices i\, . . . ,if- and ji, ■ ■ ■ ,jk at the right hand side of the 
equality (3.3) can be also written as 

k 

E E II Vfa'Vi.u) • (3-5) 

il<...<i fc ae&k r=1 
ji<---<3k 

Hence by choosing the function \P(u,v) as in (3.1) our proposition implies that 
the determinant (2.32) equals 

JV 

A( Xl ,...,x N ) (u + k)...(u + N -1) x 

k=0 

Y d\m\(xi,...,x N )m\(yi,y2,...). (3.6) 

i(\)=k 

Indeed, if \P{u,v) is the rational function (3.1) then the sum (3.5) equals 

k 



i 1 <...<i k aE&k r- 
ji<---<3k 



= 1 X i r yj<r(r) 
k 



E E E II i x ^vu {r) ) lr 



ii<...<i fe ae&k h,...,lk = l r=l 
ji<...<jk 



E E E II i^rVUr-)) 

ii<...<i fc a,re&k £(X)=k r=l 

h<—<3\ 
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Here A ranges over all partitions A of length k . The sum displayed in the last line 
equals the sum in the second line of (3.6), by using the expression (2.1) for the 
polynomial m\(xi, . . . ,xn) and a similar expression for m\(yi, 1/2, ■ ■ ■ ) ■ Thus 
our proposition implies the theorem as stated in Subsection 2.5. 



+ E V(x lt yi) 



3.2. Expanding the determinant. We will prove the proposition by induction on 
N. The case N = 1 is the induction base. Here the left hand side of (3.3) is 

00 

u 

1=1 

The right hand side of (3.3) is then the same by definition, since A(x\) = 1 . 

Now take N > 1 and assume that the identity (3.3) holds for N — 1 instead 
of N. For each index i = 1, . . . , N we will for short denote 

Ai = A(x 1 ,...,x i ,...,x N ) 

where as usual the symbol Xi indicates the omitted variable. By expanding the 
determinant at the left hand side of (3.3) in the first column and then using the 
induction assumption with u + 1 instead of u , we get the sum 

JV 



^(-ir+v-^x 



i=i 

N-l k 

(u + k + l)...(u + N-l) II *(**r,V*.) + 

fc=0 ii,...,i fe /i r=l 

31, ■■■,3k 

N 00 

^(-lr+V- 1 Yl #(*i>vi) A i >< 

i=l 1=1 
N-l k 

(u+k + l)...( U + N-l) il *( X ir>Vjr)- ( 3 - 7 ) 

fc=0 r=l 

Now consider the sum 

N 00 k 

^(-i^+V 1 ^ E E n^yi)J[n^y3r) 

i=l 1=1 ii,...,ik^i r=l 

ji, ■■■,3k 

coming from the last two lines of the display (3.7). This sum can be written as 

N k 

E(-i) i+1 ^- 1 ^, E E *(xi,yi) H #(xi r ,y3r) + 

i=l n,...,i k ^i tyji,...,j k r=l 

3i, ■■■,3k 

N k k 

E(-i) i+ V -1 ^ E E nwi.) II *(*ir,vir)- (3-8) 

i=l i\,...,i k =£i s=l r=l 
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We will prove that the sum in the second line of (3.8) equals 



N 



i+l x N-l 



^ e n 



Vjr)- 



(3.9) 



i=l 



h,...,i k ^i r=l 



Due to that equality the sum (3.7) can be rewritten as 



N 



x 



i=i 



N-l 



( E (u+k)( U +k + l)...(u + N-l) E II *( X ir,Vjr) + 



k=0 



ix-,---,ik^i r=l 

31, ■■■,3k 



AT-1 



E (u + k+ !)...(« + AT-1) £ II *0*r,V*.))- (3-10) 



ii,...,ik^i 
j lt ...,j k ^l 



r=l 



In the second line of (3.10) we can include the index k = N to the summation 
range without affecting the sum since k distinct indices i±, . . . , ik ^ i exist only if 
k < N . In the third line we can replace k + 1 with k where k = 1, . . . , N. We get 



N 



k-1 



E (u+k)...( U + N-l) E VfalVl) J[ VfarlVjr) 



k=l 



ll,...,l k -i,l 
j!,...,jk-l,l 



Here the set 

{((h,ji),---,( i k-i,jk-i),( i J)} 

can be any of the sets (3.4) appearing in (3.3), provided that in (3.4) one of the 
indices ii, . . . , ik coincides with the given index i . That one index can be then 
denoted by ik because the order of the k elements of the set (3.4) does not matter. 

These observations will show that the sum displayed in the second and the 
third lines of of (3.10) equals 



J2(u + k)...(u + N-l) E 1"! Hx^Vjr) 



k=0 



il,...,ik r=l 
ji, — ,3k 



In particular, they will show that this sum does not depend on the index i . 
Hence we will have the identity (3.3) proved, by expanding the determinant 
A(xi, . . . ,xn) in the first column. 
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3.3. Two sums. We need to prove for k = 0, . . . , N — 1 that the sum (3.9) equals 
the sum in the second line of (3.8). By using (3.2) the last mentioned sum can 
be written as 



N 



E E Ef-D^V- 1 ^ 



x 



i=l ix,...,i k ^i s—1 

jl,--;3k 



V ^ _ x II V^Vir) 



r=£s 



which in turn can be rewritten as 



N 



E E EC-^V" 1 ^ II ^y Jr ) + 



i—1 i x ,...,i k j^i s=l 
N 



r=l 



E E E(-D i+1 -."- 14 . 



x 



i=l ii,...,i fe /i s=l 

JlvvJfe 



The summands in the first line of the display (3.11) do not depend on s . 
Hence their sum equals (3.9). Let us show that the sum appearing in the second 
and the third lines of (3.11) equals zero. By opening the brackets in the third line 
and then swapping the running indices i, i s in each term coming from the second 
fraction in the brackets, the sum in the second and the third lines becomes 



N k k 

E E E(n*(*v.^))x 

i=l i x ,...,i k j^i s=l r=l 
ji,---,jk 

x"- 1 x ls Aj + (-1Y-+ 1 x^ 1 x t A la 



Here in the first line the product over r = 1, . . . , k depends neither on the index s 
nor on the choice of the index i ^ ii, . . . , ■ The fraction in the second line does 
not depend on the indices j±, . . . ,jk ■ We will show that for any fixed distinct 
indices i±, . . . , G {1, . . . , iV} the sum of these fractions over s = 1, . . . , k and 
% 7^ ii, . . . , if- is equal to zero. This will complete our proof of the identity (3.3). 

Let cr range over & n . By the definition of the Vandermonde polynomial, the 
sum of the last displayed fractions over s = 1, . . . , k and i ii, ■ ■ • ,ik equals 

* rn N-l „ N-a(i s ) + l _ N-l N-a(i s ) + l 

E E E (-i)' — x ~ x '; II 
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which is in turn equal to the sum 

k <r(i a )-l 



E E E E x?-' x*-'^- 1 n < • ( 3 - 12 ) 

i^ii,...,ife s=l <j(i) = l r=2 j^i,i s 

Recall that here N > 1 . One can easily prove by induction on iV = 2,3,... 
that the total number of terms in the sum (3.12) equals 

(N — 1)! (N — 2) (N-k)k/2. 

However, we shall not use this equality and will leave its proof to the reader. In 
the next subsection, we will prove that all terms in (3.12) cancel each other. 



3.4- Cancellations . The sum (3.12) is taken over quadruples (i, s, cr, r) . Take any 
of them such that 

cr _1 (r) ^ i 1 ,...,i k . 

Such a quadruple will be called of type I. Denote a _1 (r) = % . Observe that 
i 7^ i because o~(i) = r > 2 while a(i) = 1 . Put 

a =o-T iT (3.13) 

where G &n is the transposition of i and % . Then a(i) = r and cr(i) = 1 . 
We also have a(z s ) = a(i s ) because z s ^ i, i . Hence the quadruple (T,s, o",r) 
appears in (3.12) together with the quadruple (i,s,a, r) . But the summands 
in (3.12) corresponding to these two quadruples cancel each other. Indeed, 



v „N-r „ N-a(T) ^ N-a(i s )+r-l y[ ^N-a(j) _ 

j^i,T,i s 

o ,. A'-rrf/) A'-, /V -a(, l + r-1 ^N-a^j) 



(-l)* Xi N -*W x N-r x £-*M+ r - 1 j] 

j^i,T,i 



Note that here 

a~ l (r) = i ^ii,...,ifc 

hence s, cr, r) is also of type I . Moreover, by applying our construction to the 
latter quadruple instead of (i, s, cr, r) we get the initial quadruple (i, s, cr, r) back. 
Next take a quadruple (i, s, cr, r) showing in (3.12) such that for some index s 

a (i 8 ) - r + 1 = a(is) ■ 

Such a quadruple will be called of type II. Note that s because r^l. Put 

O = CFT i 3 i s . 

Here a(i) = a(i) = 1 because i ^ i s ,i§ ■ We also have a(i s ) = cr(i s ) . Hence the 
quadruple (i, s, a , r) appears in (3.12) together with (i, s, o~, r) . The summands 
in (3.12) corresponding to these two quadruples cancel each other. Indeed, 
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- (-lV x N ~ r x N ~ a(ia ' )+r ~ 1 x N ~ a{i "' ) TT x N ~ a{j) = 

3 ^ J ^ S J ^ S 

X - 7V " r a;. JV " <f(is) a:. JV "^ (i " )+r " 1 TT x?~°^ . 

Here 

-r + 1 = a{i s ) - r + 1 = a{i s ) 

and (i, s, a,r) is also of type II . Moreover, by applying our construction to the 
latter quadruple instead of {i, s, a, r) we get the initial quadruple (i, s, a, r) back. 

Note that the above two constructions differ for the quadruples of type 
I and II , while a quadruple can be of both types simultaneously. However, the 
summands in (3.12) corresponding to quadruples of any of the two types still 
cancel each other. Indeed, take any quadruple (i, s, a, r) of type I which also 
has type II . By applying our first constuction to it we get another quadruple 
(T,s, <r,r) of type I, where a is defined by (3.13) while % = a _1 (r) . But then 
for a certain index s we have 

a {i s ) - r + 1 = a(i s ) - r + 1 = a(i s ) 

because {i, s, a, r) also has type II . Hence the quadruple {T,s, a , r) is of type II 
too. We could similarly check that the result of applying our second construction 
to {i, s, a, r) is not only of type II but of type I as well. However, this is already 
not needed for the cancellation. So we will leave checking it to the reader. 

Finally, take any quadruple {i, s, a, r) appearing in (3.12) which is neither of 
type I nor of type II . Such a quadruple will be called of type III. Here r = a{i s ) 
for some index s because {i, s, a, r) is not of type I . For some index % ^ ii, . . . , 
we also have 

a {i s ) -r + 1 = a{i) (3.14) 

because {i, s, a, r) is not of type II . Observe that here the four indices i,ig,i s ,T 
are pairwise distinct. Indeed, here we have i,T ^ i s ,is by definition. Further, 
here i ^ % because a{i) = 1 while a{i) ^ 2 by the definition (3.14). Furthermore, 
here i s ^ i s because r < a{i s ) while r = a{i s ) ■ Put a = ar where r G &n 
cyclically permutes the indices i,i 5 ,i s ,« and leaves all the remaining indices 
fixed. More exactly, 

r: i \- >Zs(-7-z s (-7-zi— > i . 

Let f be the number at either side of equality (3.14). Consider the quadruple 
{i, s, <7, f) . Here a {%) = a{i) = 1 by definition. Due to the range of r we also 
have 

2< f ^a{i s )-l = a{is)-l. 

Therefore the quadruple {T, s, a, f) appears in (3.12) together with {i,s,a,r) . 
The summands in (3.12) corresponding to the two quadruples cancel each other. 
Indeed, we have the equality 
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- x N ~ r x N ~ a{i " ) x N ~ a{ia)+r ~ 1 x^~ a{i:) Yl x N ~ a{j) = 

(-1)* x N ~° (i) x N ~° (i ~ s)+f ~ 1 x N ~° (is) x^'^ Jl x N ~° {j) . 

Here f = a(i s ) so that (T, s,a,f) is not of type I. We also have 

a(is) - f + 1 = a(i) 

so that (z, s, a, f) is not of type II . So this quadruple is of type III. Moreover, 
by applying our third construction to this quadruple instead of (i, s, a, r) we 
get the initial quadruple (i, s, a, r) back. Thus all summands in (3.12) cancel each 
other. We have now completed the induction step in the proof of our proposition. 
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